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ABSTRACT 

In order to understand the evolution of the interstellar medium (ISM) of a galaxy, we 
have analysed the gas and dust budget of the Small Magellanic Cloud (SMC). Using the 
Spitzer Space Telescope, we measured the integrated gas mass-loss rate across asymptotic 
giant branch (AGB) stars and red supergiants (RSGs) in the SMC, and obtained a rate of 
1.4xlO~ 3 M Q yr . This is much smaller than the estimated gas ejection rate from type II su- 
pernovae (SNe) (2-4xl0~ 2 M yr -1 ). The SMC underwent a an increase in starformation rate 
in the last 12 Myrs, and consequently the galaxy has a relatively high SN rate at present. Thus, 
SNe are more important gas sources than AGB stars in the SMC. The total gas input from stel- 
lar sources into the ISM is 2-4xl0~ 2 M yr -1 . This is slightly smaller than the ISM gas con- 
sumed by starformation (~8xlO~ 2 M yr~'). Starformation in the SMC relies on a gas reser- 
voir in the ISM, but eventually the starformation rate will decline in this galaxy, unless gas 
infalls into the ISM from an external source. The dust injection rate from AGB and RSG can- 
didates is lxlO~ 5 M G yr~'. Dust injection from SNe is in the range of 0.2-1 lxl(T 4 M G yr -1 , 
although the SN contribution is rather uncertain. Stellar sources could be important for ISM 
dust (3 x 10 5 M ) in the SMC, if the dust lifetime is about 1 .4 Gyrs. We found that the presence 
of poly-aromatic hydrocarbons (PAHs) in the ISM cannot be explained entirely by carbon-rich 
AGB stars. Carbon-rich AGB stars could inject only 7xlO~ 9 M yr~' of PAHs at most, which 
could contribute up to 100 M of PAHs in the lifetime of a PAH. The estimated PAH mass 
of 18OOM in the SMC can not be explained. Additional PAH sources, or ISM reprocessing 
should be needed. 

Key words: galaxies: evolution - galaxies: individual: the Magellanic Clouds - (ISM:) dust, 
extinction - stars: AGB and post- AGB - stars:mass-loss - (stars:) supernovae: general 



1 INTRODUCTION 

Stars and the interstellar medium (ISM) of galaxies experi- 
ence a constant exchange of gas and dust. Stars are formed in 
molecular clouds in the ISM, and evolve. Elements are syn- 
thesised in stellar interior, and eventually they are ejected 
at the end of the stellar life. Gas ejected by stars is en- 
riched with metals, in comparison to with the gas where the 
stars were initially formed. Dust grains are formed around 
evolved stars and supernovae (SNe) and they are injected 
into the ISM. They might be processed in the ISM, though 
that remains uncertain. Eventually, dust grains are destroyed 
by shocks generated by SN blast winds. The lifecycle of mat- 
ter drives the evolution of the ISM, and ultimately, the evo- 
lution of galaxies. 

The concept of the lifecycle of matter, in particu- 
lar gas content, is well accepted and has been adopted 
in the chemical evolution models and stellar popula- 
tion models of galaxies (e.g. iPagel & Tautvaisienel 119981 : 
iBruzual & Chariot! 120031; iKodama & Arimotolll997h . How- 
ever, it had been difficult to actually measure the gas and dust 
feedback from stars into the ISM in real terms. The Spitzer 



Space Telescope dWerner et al.l2004h has provided an oppor- 
tunity to meas ure the gas feedback from stars in the Mag- 
ellanic Clou ds dMatsuura et alJl2009llSrinivasan et alj|2009t 
iBover et alj|2012l) . delivering new constraints on chemical 
evolution models. 

Dust is one of the important contents of galaxies. Dust 
absorbs energy emitted from stars within a galaxy, and re- 
emits it in the infrared to sub-mm wavelengths. The pres- 
ence of dust controls the energy input and output from the 
ISM of galaxies, and affects the spectra l shape of the galaxy 
and the underlinin g physics of the ISM (Gall iano et al.. 2005: 
iDunne et alj|201 ll) . 

Despite such an important role in the physics of the 
ISM and galaxies, it is sti ll not well establi s hed how dust 
mass in the ISM evolves dSloan et alj 120091: iDraind 120091 
iTielens et all 120051: ICalura et all 120081) . This requires de- 
tailed studies of the dust formation and destruction pro- 
cesses in stars and the ISM. Dust grains are considered 
to be formed in a wide range of object s, part i cularl y 
stars in the late phase of evolution (e.g. iGehrzj 1 1 9891). 
Such theoretical studies have been conducted (Nozawa et al. 
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20031; iMorgan & Edmunds! 120031: iFerrarotti & Gaill 120061: 
Zhukovska et al. l2008t IValiante et al.ll2009l) but are poorly 



constrained by observations at present. Once grains are 
ejected from stars into the ISM, it has been proposed that 
dust grains are re-processed and destroyed. Dust grains 
might grow in the ISM, using stellar dust as seeds (e.g. 
iTielens et alJl2005t lDraind l2009). but it is more difficult to 
directly measure these processes. This paper focuses on es- 
tablishing an understanding of what the contribution of stel- 
lar dust is to the evolution of dust in the ISM. 

One of the important dust sources are asymptotic giant 
branch (AGB) stars, low- and intermediate-mass (1-8 M ) 
evolved stars. It has been well esta blished that du st grains are 
formed in the AGB outflow (e.g. lHabind ll996). However, 
it is challenging to make a quantitative analysis; namely to 
measure dust formed in evolved stars in entire populations 
across galaxies, and evaluate their contribution to dust in the 
ISM, because this requires a large survey of AGB stars in 
a galaxy. One of the pioneering studies was based on the 
IRAS survey of the solar neighbourhood dJura & Kleinmannl 
1989), and the Spitzer Space Telescope open ed up a possibil- 
ity for such studies beyond the Milky Way ( Ma tsuura et alj 
120091: ISrinivasan et alj|2009h . 

High-mass stars are considered to form dust in vari- 
ous evolutionary stages, such as red supergiants (RSGs), 
Wolf-Rayet stars, luminous blue variables (LBVs), and 
supernova (SNe). Due to limited number of sample, the 
mass dust formed in high mass stars, and in particu- 
lar, SNe, is uncertain. Recent Spitzer and Herschel stud- 
ies of SNe and SN remnants (SNRs) show that type-II 
SNe can form dust, but the reported dust masses range 



from 10~ 4 to ~lMn dSugerman et al.1 l2006t iMeikle et al 



2007|;lBarlow et al.l2010HMatsuura et alj201 Ubomez et al 
2012). Type-la SNe, which have an origin in low- and 



intermediate-ma ss binary stars, app ear to form very small 
amounts of dust (Gomez et al. 2012b- 

The explosion of high mass stars creates expanding 
winds, which of ten have multiple velocity components 
(Kia?r ^t~ai1l2010l) . When the fast blast wind collides with 
the ISM dust, ISM dust grain s can be destroy ed via sputter- 
ing and shattering processes dJones et al.ll 1 996b . The lifetime 
of dust is determined by its destruction by SN shocks. As- 
suming homogeneous gas and dust distributi ons, the lifetim e 
of dust is estimated to be nearly 1 Gyr dJones et al.| [l996). 
Gas and dust are distributed inhomogenously in the ISM, so 
as it is very difficult to estimate a precise dust lifetime. 

In this paper, we measure the gas and dust injection rate 
from AGB stars and RSGs in the Small Magellanic Cloud 
(SMC ). This galaxy is only 56kpc away dMcCumber et al.l 
2005) and is a molecular-poor galaxy, yet has some ongo- 
ing starformation. The measured gas and dust masses ejected 
from AGB stars and RSGs are compared with masses ejected 
from SNe, and with the starformation rate in the SMC. We 
discuss the current problems in understanding the gas and 
dust budget, and its evolution in these galaxies. 



2 MASS-LOSS RATES OF AGB STARS AND RED 
SUPERGIANTS 

Mass-loss rates of AGB stars and red supergiants correlate 
well with infrared excess. Often near-infrared photometric 



points are used to represent the photospheric flux, while 
mid-infrared fluxes illustrate the emission from the dust 
thermal emission. Taking near- and mid-infrared colours is 
an indicator of mass- loss rates dLe Bertre & Winters! 1 19981; 
IWhitelock etaD 19941) . 

The correlation between colour and mass-loss rates 
have been established for LM C and SMC carbon-rich AGB 
stars (Matsuura et al.l (12009). who used the mass-loss mea- 



surement s fromlGroe newegen et al. (20 Q7|), and pho t ometr ic 
data from lMeixner et al.ld2006l) and lSkr utskie etaD {2006)). 
The resultant relations were 

log M g = -6.20/(([3.6] - [8.0]) + 0.83) - 3.39 (1) 

in the range of 1 < [3.6] - [8.0] < 9 

log M g = -14.50/(CSr s - [8.0]) + 3.86) - 3.62 (2) 

in the range of 1 < K s - [8.0] < 9, where M g is the gas 
mass-loss rate. 

The relation for K s - [24] for the same data set is found 
to be 



logMg = -20.25 l((K s - [24]) + 4.98) - 3.47, 



(3) 



We further derive the r elatio nship for oxygen-rich stars, 
using Groenewegen et al.l (l2009l) 's analysis of oxygen-rich 
AGB stars and red supergiants in the Magellanic Clouds. 
We c orrelate their photom etric data with Spitzer and 2MASS 
data (Gordon e t al.ll20lTl) . The resultant correlation is found 
in Fig[T] and the fits to these data are given as 

log M g = -18.97/(([3.6] - [8.0]) + 2.698) - 0.9954 (4) 

in the range of 1 < [3.6] - [8.0] < 3 

log Mg = -U.2S/((K S - [8.0]) + 2.297) - 2.888 (5) 

in the range of 1 < K s - [8.0] < 7, and 

log Mg = -38.96/((K s - [24]) + 4.903) - 1 .522, (6) 

in the range of 1 < K s - [24] < 9. 

It is difficult to measure the mass-loss rates of stars with 
very little infrared excess and blue colour cut-offs. How- 
ever, contributions from very low-mass loss rate stars are 
not important for the overall global gas a nd dust budget 
dLe Bertre et al.ll2001l;lMatsuura et a l. 2009). 



3 OBJECT CLASSfflCATIONS 

The first step is to understand the infrared photometric sur- 
vey and classify the point sources. Colour magnitude di- 
agrams (CMDs) and colour-colour diagrams ( CCDs) are 
commonly used to classify the point sourc es (e.g. lBlum et al.l 
120061: llta et all 120081: IWoods et al.ll2011l) . In this work, we 
use the LMC objects to set object classifications based on 
a CMD and a CCD, and apply the classification method to 
SMC objects. 

In our previous study, we focused on classification 
based on [3 .6]-r8.01 v.s. [8.0] CM D to extract carbon-rich 
AGB stars (Matsu ura et al.l l2009t) . We are going to revise 
this classification, and particularly optimise the classifica- 
tions for high mass-loss rate AGB stars, which are dominant 
for the gas and the dust budget. 
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3.1 Cross-identifications 
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We start the object classification process by cross- 
identifying the LMC photometric data with spectroscopi- 
call y known objects in the LMC. 

Mats uura et al.l (I2009I) have assembled the LMC objects 
with spectroscopical classifications and we use these sam- 
ple: Sample of M-type giants and supergian ts, carbon-rich 
stars, p l anetary nebu l ae are as sembled from | Kontizas et al.l 
(2001 b. ICioni etal l d200ll). ISanduleak & Philid dl977l). 
iBlanco et all d 19801) i IWesteriund et al l d 19811) . IWoodetal. 
(119831 ) . IWood et al.l d 19851) . iHughesI d 19891) and iReid et al. 
(1988 ). We added sam ples of Wolf Rayet (WR) stars 
(Brevsach er et alj 1 19991) and S Dor variables (Lumi 

nous Blue Variabl es (LBVs) ; |Van Genderen 2001 

post-AGB stars dGielen et al.l 120091: IVolk etalJ 1201 
Matsuura et alj|201 lb. and planetary nebulae (Wood s et al.l 
201lHMatsuuraetalJl2011l) . 
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Figure 1. Correlation between mass-loss rate (in M e yr~') and colour for 
oxygen-rich AGB stars and RSGs 



3.2 Object classifications 

3.2.1 LMC colour-magnitude diagram 

Figure|2] shows the [3.8]-[8.0] vs [8.0] CMDs of the LMC 
objects. This diagram was used for object classifications of 
carbon-rich AGB stars and oxy gen-rich AGB stars/re d su- 
pergiants in our previous study dMatsuura et al.l l2009). The 
black lines shows the separation of th ese two types . LMC 
Spitzer spectroscopic observations ( Kem per" et al.l 1 2010) 
confir m that this classification is effective dWoods et al.l 
l20Tlh . 

This CMD classification is simple and largely correct. 
Actually, contamination of a small number of red objects 
(carbon-rich post-AGB stars and PNe and distant galaxies) 
into AGB stars and RSGs could potentially change the anal- 
ysis of the global gas and dust budget, which is the final 
aim of this paper. The further and more severe problem is 
contamination of a few, but high mass-loss rate oxygen- 
rich AGB stars which fall into the carbon-rich AGB region 
at about [3.6]-[8.0]~2 and [8.0]~6mag. We introduce one 
more step in the analysis, as described in Sect. 13.2.21 to min- 
imise these contaminations. 

Many point sources are found in the region between 2 < 
[3.6] - [8.0] < 4 and [8.0] < 12. This is ass ociated with 
distant galaxies (Kozlowski & Kochanek 2009). 



3.2.2 LMC colour-colour diagram 

In order to improve object classifications further, we intro- 
duce another colour-colour diagram, the K-[8.0] vs K-[24] 
of the LMC objects, as plotted in Fig . [3] . 

Prior to the plot, point sources as sociated with distant 
galaxies dKozlowski & Kochanekl l2009) were removed. We 
defined their distribution as [3.6] - [8.0] > 2.0 and [8.0] < 
10 from Fig. [2] These galaxies would appear approximately 
between 1.5 < K - [8.0] < 2.5 and 5 < K - [24] < 7 in the 
CCD, if they had been plotted. 

We overlay the mass-loss vs colour relationship derived 
in Section top of the CMD. There are two separate se- 
quences for oxygen-rich and carbon-rich objects. The red- 
der the colour, the higher the mass-loss rate, in general, as 
the equations show. 
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Figur e 2. The infrared colour-magnitude diagram of the LMC objects. Dots represent all point sources obtained from the SAGE catalogue (Meix ner et alj 
2006). Circles show the spectroscopically identified carbon-rich AGB stars. Upword-pointing triangles show oxygen-rich AGB stars and red supergiants 
(RSGs). In addition to the AGB stars, S Dor type variables (luminous blue variables), and Wolf Rayet (WR) stars are plotted. The dashed lines show the 
expected detection limit of the final SAGE catalogue. 



We use this CCD and the CMD to classify objects. We 
cross check the reliability of these classifications with SIM- 
BAD. We found that the contaminations are not a big issue 
for the LMC gas and dust budget, but it does change slightly 
(10%) for the SMC. A cross-check with SIMBAD found that 
non-AGB red objects that contaminated the AGB star s are, 
if they have kno wn identifications, P Ne dHenizdl 1955b and 
YSO candidates dBolatto et alj|2007l) . 



stars in the LMC. We take into account the difference in 
the number of stars in these two galaxies. Taking the ab- 
solute V-band magnitude as a measure of number of stars 
in a galaxy, the SMC should have about 7 times fewer stars 
than the LMC dWesterlundll 19901) . The difference in the red 
carbon-rich AGB stars appears be slightly larger than the dif- 
ference in the number in the stars, but this is not conclusive. 



3.2.3 SMC colour-magnitude diagram 

Once we establish the object classifications using the LMC 
objects, we apply the same method to the SMC objects. 

Figure|4] shows the CMD of SMC poi nt-sources ex- 
tracte d from the SAGE-SMC catalogue (iGordon et alj 
1201 lb . The black lines show the separation between carbon- 
rich AGB stars and oxygen-rich AGB stars/red supergiants, 
which was derived from LMC objects, but are scaled by 
the difference of distance moduli between these two galax- 
ies. The distance moduli of t he LMC and the SMC are 
18.5 and 18.9 m ag respectively (Nikola ev & Wein berg 2000; 
Westerlund 1990), so that the separation line was shifted by 
0.4 magnitude fainter in [8.0] for the SMC objects. 

There is a clear difference in LMC and SMC CMDs: 
fewer red stars are found in the SMC than the LMC. Simi- 
larly, there are only three carbon-rich candidates with very 
red colours ([3.6]-[8.0]>5.0), whereas there are 33 such 



3.2.4 SMC colour-colour diagram 

Figure[5] shows the same combination of CCD as Fig. [3] but 
for the Small Magellanic Cloud (SMC). There is a slight dif- 
ference in the distributions of stars compared with the LMC 
diagram: the SMC has very few oxygen-rich AGB stars with 
colour redder than K - [8.0] > 1.8. This corresponds to ap- 
proximately a mass-loss rate higher than 10~ 6 M yr . 

The SMC CCD shows that there are a reasonable num- 
ber of stars found in region (a), which are carbon-rich AGB 
candidates. Oxygen-rich AGB stars, represented by region 
(b), are relatively scarce. There are even fewer oxygen-rich 
AGB stars in the SMC than the LMC (Fig|3). 
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Figure 3. Colour-colour diagram of the point sources in the LMC. The objects with spectroscopic classifications are marked with symbols. We divide the 
colour-colour diagram into 4 regions (a)-(d) and dominant types of objects in these regions are indicated. Oxygen-rich stars and carbon-rich stars follow two 
separat sequences. Two dashed lines show the estimated mass-loss rate and colour relations for oxygen-rich and carbon-rich stars, respectively, and they are 
calculated from equations (l)-(6). The numbers alongside the lines show the mass-loss rates in log-scale (Sect|2). 



Table 1. Gas and dust injected into the ISM of the LMC and the SMC 
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SMC 
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Gas 


Dust 


Gas 


Dust 




(10- 2 M Q yr-') (10" 


"SMoyr- 1 ) 


(lO^Moyr- 1 ) (10" 


^Moyr- 1 ) 


Carbon-rich AGB stars 


0.7 


4 


0.08 


0.4 


Oxygen-rich AGB + RSGs 


0.8 


4 


0.06 


0.3 


Type II SNe 


6-13 


7-400 


2-4 


2-110 


WR stars 


~ 0.1 




~ 0.01 




OB stars 


0.1-1? 




-0.03-0.3 




Star-formation rate 


20-30 
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4 GAS AND DUST INJECTION RATE FROM EVOLVED 
STARS 

4.1 Integrated mass-loss rate from AGB stars and red 
supergiants 

Combining the information from object classifications and 
mass-loss rates from the previous sections, we can estimate 
the integrated gas and dust mass-loss rates from AGB stars 
and red supergiants. This is the total gas and dust mass in- 
jected from evolved stars into the ISM. 

In order to estimate mass-loss rates of individual AGB 
stars and RSGs, we use the correlation of [3.6]-[8.0] with 
mass-loss rate. An alternative choice is K s — [8.0], but some 



high mass-loss rate stars are not detected at 2MASS K s - 
band, and these stars could contribute a significant fraction 
to the integrated mass-loss rate from evolved stars. 

Object classification is based on [8.0] vs. [3.6]-[8.0] 
colour, and we supplementary use /T,-[8.0] vs ^.,-[24] to 
remove the contamination of distant galaxies in the red 
carbon-rich AGB region, as well as high mass-loss rate 
oxygen-rich AGB stars and RSGs in the carbon-rich AGB 
region. 

By integrating the mass-loss rates of individual AGB 
stars, we can estimate the total gas and dust inputs from 
AGB stars and red supergiants into the ISM of the SMC, 
as summarised in TableQ] The contribution from oxygen- 
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Figure 4. The colour-magnitude diagram of the SMC objects. The regions which are expected to be dominated by carbon-rich AGB stars and oxygen-rich 
AGB stars and red supergiants are indicated. 



rich AGB stars and RSGs contribution is 0.06x 10 -2 M yr~' , 
while carbon-rich AGB stars contribute 0.08xl0 _2 M G yr _1 
in gas. The total gas injection rate from AGB stars and red 
supergiants is O.lxlO _2 M yr _1 for the SMC. 

The estimated total gas input from evolved stars is 
2xl0" 2 M o yr _1 for LMC. This includes both oxygen-rich 
and carbon-r ich AGB stars and RS Gs. In our previous work 
on the LMC (Matsu ura et al.ll2009T) , the integrated mass-loss 
rate from oxygen-rich AGB stars and RSGS was a crude es- 
timate based on previously known surveys before the SAGE. 
Tthe current estimation has a slightly higher value. The 
integrated mass-loss rate is added up, using the contribu- 
tion from carbon-rich AGB stars (0.7xl0~ 2 M o yr _1 ) and 
oxygen-rich AGB stars and RSGs (O.8xl(T 2 M yr -1 ). 

In the SMC, dust inputs, which are directly measured 
from mid-infrared observations, are 0.4xlO- 5 Moyr _1 for 
carbon-rich AGB and 0.3xl(T 5 M yr -1 for oxygen-rich 
AGB and RSGs. 

The gas-to-dust mass ratio is assumed to be 200 for 
all A GB stars and RSGs as in iGroenewegen et all d2007l 
2009). This is the largest uncertainty in our estimate of the 
integrated gas mass-loss rate. The ratio of 200 is a widely 
assumed value for the Galactic AGB stars and RSGs, but 
it is not clear that this ratio is applicable to extra-Galactic 
evolved stars, where the metallicities of galaxies are differ- 
ent. The gas-to-dust ratio in the ISM and evolved stars is left 
for future investigations, such as ALMA observations. 

iBover etail d2012l) has estimated that the total dust in- 
jection rate is 1 x 10~ 7 M yr _I from carbon-rich AGB stars, 



6x 10~ 7 M yr~' from "extreme-" AGB stars, which they ex- 
pected to be mostly carbon-rich, and 0.7 x 10~ 8 M yr~' from 
oxygen-rich AGB stars, and 3xlO~ 8 M yr -1 from RSGs 
We re moved the ambiguity of "extreme" from iBover et"ai1 
d2012l) 's classification. We found that the contribution from 
oxygen-rich AGB st ars and red s upergi ants are much larger 
than their estimate. IWood et all d 19921) found three OH/IR 
stars in the SMC. They have not estimated the mass- 
loss rate of these particular three objects, but the typi- 
cal range of mas s-loss rate of OH /I R stars are 3 x 10 ~ 5 - 
3 x 1Q- 4 vr^ dWood et alJll992l) . I Wood etaLld 19921) es- 
timated these mass-loss rates based on infrared observa- 
tions, and they adopted a gas-to-dust ratio of 200. The 
sum of the dust mass-loss rates from these three OH/IR 
stars alone is estimated to be 2 x 10~ 7 -2 x 10 _6 M Q yr _1 . 
iGroenewegen et ail d2009l) analysed the mass-loss rates of 
AGB stars and RSGs in Magellanic Clouds, and 1 1 MSX 
named oxygen-rich SMC stars have a total dust mass loss 
rate of 4 X lO _7 M yr _1 . These values show that our es- 
timate s are reasonable. This does not affect IBover etail 
d2012l) 's overall conclusion that AGB stars and RSGs are 
not the dominant source of dust in the ISM, which we also 
found. 



4.2 SN gas and dust injection rate 

In order to estimate the average gas and dust injection rate 
from core-collapse SNe into the ISM, we need to estimate 
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SN rates and the progenitor mass, as well as dust mass 
formed in the SNe themselves. 

The SN rates have been estimated for both the LMC and 
the SMC, using the number and age of the supernova rem- 
nants (SNRs) in thes e galaxies as a starti ng point. In our pre- 
vious LMC_study (Matsu ura et alj|2009l). we used SN rate s 
fro m iMathewson et all d 19831) and iFilipovic et al.l dl998l) . 
The IMathewson et al.l dl983D est i mate is a facto r of tw o 
lower than IFilipovic etalj d 19981) . IFilipovic etail dl998l) 's 
study is more recent, so we adopt th is value here. Their rate 
includes both type la and II SNe. Tsuiimoto et"ai] d 19951) 
estimated the ratio of type la over type II is 0.2-0.3 for 
both the LMC and th e SMC. Combining the SN rate from 
Filipovic et al.l (fi~998|) and type la and type II ratio from 



Tsuiimoto et al.l d 19951) . we obtain type SN II rate of one 



in every 125-143 years in the LMC and one in every 438- 

500 years in the SMC. 

In our previous study of the LMC (Matsuura et al. 

2009), we have used the average SN progenitor mass of 8 M 
from SN progenitor surveys in the galaxies dSmartt et al.l 
120091) . Eight solar masses correspond to the lowest range of 
mass that ends their stellar evolution as SNe. One possibility 
is that this average progenitor mass might have some bias 
towards the lowest end of SN progenitors, because they are 
more numerous than higher mass stars and more frequently 
detected. We take this possibility into consideration. 

A mean mass of high-mass stars is much higher than 
8M n ; if we ta ke a mean mass of the initial mass function 
dKroupall200il) with the high mass cut off of 50 M , the av- 



erage mass of high-mass stars would be about 16 M . We 
take this to give a maximum possible mass of gas ejection 
rate per SN. 

Considering the progenitor mass range and SN rate, the 
integrated gas inputs from type II SNe are approximately, 
6-13xlO- 2 M Q yr _1 for the LMC, and 2^xl0- 2 M o yr 1 for 
the SMC. 

SNe form dust in their ejecta, however, measured 
dust masses from SNe and SNRs have a wide range. 
Mid-infrared observations found the lower limit of dust 
mass to be about ~ 10~ 5 M^ dMeikle et al.l 120071). bu t 
up to ~ 10~ 2 Mm dSugerman et al.l 120061: iFox et all l20lTh . 
Recently, the Hersc hel Magellanic Survey (HERITAGE; 
iMeixner et all [2 010) found a large amount of dust in 
SN 1987A dMatsuura et al.ll201 ll) . and its mass is about 0.3- 
0.7 M , showing SNe can form significant dust mass, fol- 
lowing a measurements in Galactic SNR, CasA (0.08 M 
of dust; iBarlow et al1l2010HSibthorpe et al.ll2010l) . It seems 
that the measur ed dust mass does not correlate with the pro- 
genitor mass dOtsuka et alj 1201 It iGall et al]|201 lb . In our 
estimate of the SN dust injection rate, we take two cases, 
~ 1O~ 2 M as a conservative c ase, as taken in our previ- 
ous study (Matsuur a et al J 120091) . The other case is 0.5 M 
per SNe. This assumes that the dust detected in SN 1987 A 
after 25 years from the explosion would not be destroyed 
in the later stage by the ISM and SN wind collisions. The 
dust mass per event brings the dust ejection rate to 7- 
400x10 " 5 M yr' for the LMC, and 2-llOxlO^Moyr 1 
for the SMC. 
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4.3 Other dust and gas sources 

There are other possible sources which could contribute to 
the gas and dust budgets, such as Wolf-Rayet (WR) stars, 
luminous blue variables (LBVs), novae and OB stars. We 
briefly discuss their contributions to the budget. 

Studies of Galactic WR stars show that late-type (WR9 
and W C 10) stars form carbonaceous dust dWilliams et al.l 
Il987l) . In the Magellanic Clouds, such late-type WR stars 
have not been found yet. Among known WR stars, the LMC 
WR star (H P 36402; WC 4 ) was detected at mid-infrared 
wavelength (IWilliamsll20i~il) . and its mass is estimated to be 
1.5 x 10~ 7 M Q . It is not clear what the time scale of this dust 
formation and ejection rate is, but WR dust input is much 
smaller than that of AGB stars and SNe combined in the 
LMC. There is no report of dust formation in SMC WR, so 
far. 

There are just ov er ten WR stars known in the SMC 
(Paseman n et al.l [201 ll) . Their mass loss is line d riven, and 
the mass-loss rate decreases w ith metallicity dKudritzkil 
12002b ICrowther & Hadfieldll2006l) . Assuming eac h WR star 
ejects approximately 1 x lO~ 5 M yr~' of gas dCrowfheri 



1 per star (Prinia 1987), the total mass-loss rate 
3 x 10~ 5 -~ 3 x lO^Moyr 1 . In the LMC, 



2000), WR stars eject is about 1 x 10~ 4 M yr~' of gas into 
the ISM in the SMC. T he LMC WR catalog lists 134 stars 
dBreysacher et alJI 19991) and about 1 x lO~ 3 M yr~' of gas 

is ejected from WR stars into the ISM in the LMC 

There are over 300 OB stars in the SMC dLamb et all 
1201 lb . If their mass -loss rate is typically about 10~ 8 
10- 6 M G yr ■' 
is about ~ 

the total mass-loss rate from approximately 1000 OB stars 
dBastian et al.ll2009h is about ~ 1 x 10 4 -~ 1 x lO^Moyr -1 . 

Luminous blue variables (LBVs) coul d potentially 
make significant amount s of silicate dust dMorris et al.l 
Il999t iGomez et alj 120101) . Among the SDor star catalog 
of I Van Gendered d200ll) . three LBVs have an infrared ex- 
cess in the LMC, w hich are R71, R127 and BAT 99-83 
(Bonanos et al. 200 9)), indicating dust formation in these 
stars. In the SMC iBonanos et all d2010h found three LBVs 
that indicate an infrared excess. However, we found that 
solely IR data can not distinguish between RSGs and LBVs 
(FigO, so that LBVs are combined together with RSGs in 
the overall gas and dust budget. 

Classical novae can form dust and observations of 
Galactic novae reported dust mass ranges between 10 ~ 10 - 
10~ 5 M Q depending on the objects dGehrz et alJll998l) . As 
far as we are aware, there is no report of dust formation in 
novae in the Magellanic Clouds, so that we do not include 
novae in the dust budget, at this moment. 

iKastner et alJ d2006l) found two LMC Be stars to have 
dust in discs. These are probably formed during occasional 
outbursts. The estimated mass-loss rate was 5 x 1O~ 4 M yr~' 
giving the minimum mas s injection rate from B e stars into 
the ISM of 10~ 3 Mm yr~' . IBonanos et all d2010l) have cross 
correl ated optical spectr al surveys of massive stars in the 
SMC dEvans et alj|2004l) . and found that five B[e] stars have 
infrared excess. They could have contributed approximately 
1 x 10~ 5 M o yr _1 of dust. 

Two R CrB stars have an in frared excess. One is MSX 
LMC 439 dTisserand et al.ll2009h and the other is MS X LMC 
1795 dSoszvnski et al.ll2009l) . IClavton et alj d2011l) studied 
another LMC R CrB star, H V 267 1 . In our estimate, we in- 



tegrate R CrB stars into carbon-rich AGB stars, as they can 
not be distinguished from IR data only (Figf5|. 

AGB stars have one more chemical type, on top 
of oxygen-rich (M-type) and carbon-rich (C-type): S-type 
stars. They have a carbon-to-oxygen ratio almost equal to 
unity. Their dust composition is diverse: FeO, FeS, silicate 
dZiilstra et alj2004tlSmolders et alj2012b . The number of S- 
type stars is small (~1 % in the Galaxy) among AGB stars, 
so their contributions to the total gas and dust budget is neg- 
ligible. 



5 DISCUSSIONS 

5.1 Global gas and dust budget of the SMC 

In the SMC, the total gas injection rate from AGB stars, 
RSGs and SNe is 0.02-0.04 M yr _1 . This largely relies on 
the gas ejected from SNe. 

We compare this gas injected from SNe and evolved 
stars into the ISM with gas consumed by star forma- 
tion. [KennicuiljfcHodil {l986|) estimated the current SMC 
starformation rate of O.O8M yr _1 . In the SMC ISM, gas 
consumed by the starformation exceeds the gas injected 
from evolved stars and SNe (0.02-0.04 M yr -1 ) into the 
ISM. This is similar to what has been found in the LMC, 
though the deficit in the LMC is much larger (about 0.1- 
O.2M0VT- 1 ). 

The majority of the ISM gas in the SMC is present in 
the form of Hi and its ma ss is 4 x 10 s M^ dSta nimirovi cTet al.l 
ll999tlBolatto et alJ[20TH) . A large ISM reservoir could sus- 
tain the star formation, despite a deficit in the gas injection 
rate, at least on a few Gyr time scale. Unless there is a gas 
infall from an external source, the starformation rate will de- 
cline. 

Similarly, SNe and AGB stars are important dust 
sources in the SMC. SNe could contribute to a dust injec- 
tion of 2-1 lOxlO 5 M yr~' into the ISM, while AGB stars 
and RSGs contribute about lxl0~ 5 M yr~' of dust mass. 

In the Milky Way, the lifetime of dust is estimated to be 

6 X 10 8 years and 4 x 10 s years for carbonac eous and silicate 
dust grains, respectively dJones et aT]|1996l) . The lifetime of 
dust grains is constrained by their destruction by fast SN 
shocks and the d ust lifetime co rrelates with the SN Ia+II rate 
per surface area dDwekH l998). The Galactic SN rate i s esti- 
mated to be one in every 50 years dDiehl et al.l l2006). This 
rate includes type Ib/c and II, but in practice type II dom- 
inates. The type la rate in the Milky Way is about 0.3 per 
century dMatteucci et alj|2009l) . so that the total SN la and II 
rate is about one event in every 45 years in the Milky Way. 
The SMC has a rate of one event every 350 y ears, and in the 
LMC the rate is one event every 100 years (Filipov ic et al.l 
1998). To calculate the SN rate per surface area, we use the 
half-light radius (ri/2) of a ga laxy as a measure of the sur- 
face area, which we take from lTolstov et al.l d2009l) . The es- 
timated dust lifetime would be 5 x 10 8 years and 3 x 10 8 years 
for carbonaceous and silicate dust grains, respectively in the 
LMC, and 17 x 10 8 years and 11 x 10 8 yeas in the SMC. 
In the SMC, the overall lifetime of the dust would be about 
1.4Gyrs. This assumes an almost constant SN rate and star- 
formation rate over such a long period, which will be dis- 
cussed in Sect.|521 In the SMC, about 4 x 10 4 -1 x 1O 6 M 
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of dust has been accumulated over the 1 .4 Gyr history of the 
SMC. 

The estimated dust mass in the SMC ISM is 3 x 10 5 M Q 
dLeroy et al.ll2007l) . If the dust lifetime in the SMC ISM is as 
long as 1 .4 Gyrs on average, dust present in the SMC ISM 
may mainly originate from stellar dust. 

The obvious uncertainty in this estimate is the lifetime 
of dust. The half-light radius is an indicator of the galaxy 
size, but the LMC and the SMC are irregular galaxies, and 
the actual SN rate per volume would not be so simple to esti- 
mate. Furthermore, SNe might affect dust destruction in the 
SN vicinity only, and dust in the remaining regions might 
survive. The uneven distribution of dust could bring diverse 
life-time of dust. It seems that the dust lifetime largely de- 
pends on the local condition of the galaxies, and more so- 
phisticated calculations are needed. 

It has been proposed that dust condensation in molec- 
ular clouds could be a dominant source o f ISM dust in 
the Milky Way or high-red shift galaxies (iTielensi 120061 : 
Draine 2009; Mattsson 2011). These galaxies probably con- 
tain many molecular clouds within. The SMC is known to 
have a low molecular content ({israelii 1 997h . so that dust con- 
densation in the molecular clouds may not be as high as in 
the Milky Way. 

The key to understand the evolution of ISM dust is 
to understand the formation and destruction processes im- 
posed by SNe. It is now pausible to determine dust formation 
and destruction in SNe, using the current missions/projects, 
such as Herschel, SOFIA and ALMA, but also fu ture space 
missions (JWST and SPICA; iTanaka et al]l2012l) . Also we 
need to find constraints for grain growth in molecular clouds 
based on observations. 



5.2 Comparison between the LMC and the SMC 

In the SMC SNe are a more important gas source to the ISM 
than AGB stars and RSGs. The difference between them is a 
factor of 14-29. In the LMC SNe have higher gas feedback 
rates than AGB stars and RSGs but the difference is only by 
a factor of 4-9. AGB m ass-loss rates have uncert ainties of 
a factor of three at least dGroenewegen et aT]|2007l) . and the 
gas-to-dust mass ratio could have a factor of three uncertain- 

I II ,A 

ties ( Ivan Lo on 2000). The LMC difference is negligible, but 
the one in the SMC is not. There is a difference in the main 
gas sources in the LMC and SMC. The SMC has more gas 
feedback from high-mass stars th an the LMC. 

The initial mass function dKroupal 1200 lb shows that 
equivalent mass should be distributed to high-mass stars 
(>8 M ) and low- and intermediate-mass stars (1-8 M Q ). At 
the event of stellar death, high-mass stars could return most 
of their mass into the ISM. Many low- and intermediate- 
mass stars end th eir lives as white d warfs, and their masses 
are about 0.6 M jWeidemannll 19901) . thus 40-92 % of mass 
is returned to the ISM. During a constant starformation rate 
period, the ratio of ISM gas feedback from high-mass stars 
against low- and intermediate-mass stars should be about 
3:2. The measured ratio shows that SMC has an excess of 
gas feedback from high-mass stars beyond the IMF. 

The LMC and SMC ha ve experienced more or less sim - 
ilar starformation histories (lHarris & ZaritskvllIOOl [2009), 
but the SMC has an enhanced starformation history by a fac- 



tor of four in recent times in the last 12Myrs. That could 
result in a more efficient gas injection rate from SNe. 

Additionally, the lower metallicity can cause lower dust 
driven wind s, at least for oxygen-rich AGB stars and re d 
supergiants dBowen & Willsonl 1 1 99H iMarshall et~aT]|2004l) . 
The metallicities of the LMC and the SMC are about half 
and the quarter of Solar dMonket alj|1988l) . That could be 
another reason that there is a more efficient SN gas feedback 
in the SMC, because of the relatively smaller contribution of 
oxygen-rich AGB stars and red supergiants. 

5.3 PAHs in the SMC 

Carbon-rich AGB stars are con sidered to be the source 
of PAHs found the ISM dAllamandola et alJ 119891: 
ICherchneffetail Q992). PAHs are formed using C2H2 
as a parent molecule in chemical reactions. In the circum- 
stellar envelope of carbon-rich AGB stars, carbon atoms are 
tied up in CO first, and use up all oxygen, and the excess 
carbon formed carbonaceous dust and carbon-bearing 
molecules. 

From our measured carbon-rich AGB gas injection rate, 
we can estimate the upper limit of PAH mass injected from 
carbon stars. The gas injection rate from carbon-rich AGB 
stars is O.7xlO _2 M yr _1 . The C2H2 fracti onal abundance 
is app roximately 10~ 5 in the Milky Way dHowe & M illar 
119961) . This can be has a factor of few higher, depending 
on how the metallicities of the galaxies affect the amount of 
excess carbon in the stars (Matsuura et al. 2005), and to be 
predicted to be about 10~ 4 dWoods et alj[2012l) in the Magel- 
lanic Clouds. The maximum PAH injection rate, if all C2H2 
is eventually converted into PAHs, is O.7xlO _6 M yr _1 . 

The average fractional abundance of P AHs in the SMC 
is 0.6 % with respect to the total dust mass (ISandstrom et al.l 
120101) . The estimated PAH mass in the SMC to be about 
18OOM . 

The lifetime of PAHs is shorter than amorphous carbon. 
In the Milky Way, the estimated life time of PAHs is 1.4- 
1.6xl0 8 years dMicelotta et al.ll2010l) . In the SMC, the SN 
rate is lower, and the lifetime of PAHs can be as long as 
a few 10 8 years. The PAH injection rate from carbon-rich 
AGB stars is so low that the expected PAH mass from AGB 
stars would be much lower than 100 M . PAHs in the SMC 
require efficient in the ISM or in-falling into the SMC. PAH 
formation needs ISM processing. 



6 CONCLUSIONS 

We have measured the total gas and dust injection rates from 
AGB stars, red supergiants, and also estimated these rates 
from supernovae. The total gas injection from stellar deaths 
is about 2-4xlO _2 M yr _1 into the ISM. This is slightly 
smaller than the current gas consumption in the ISM by 
starformation, which is ~8xlO~ 2 M yr~'. The galaxy has a 
large gas reservoir the moment, so that it can sustain high 
starformation rate at present. Eventually the starformation 
rate is going to decline, unless an external source provides 
gas infall into the SMC. 

AGB stars and red supergiants are important sources of 
dust in the SMC. Dust production in SNe is largely uncer- 
tain, so is the lifetime of dust. Within the current uncertain- 
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ties in quantities, dust present in the ISM can be explained 
as being stellar in origin. 

PAHs in the ISM can not be explained to have originated 
only from carbon-rich AGB stars. The lifetime of PAHs is 
too short, compared with the supply from carbon-rich AGB 
stars. PAHs require formation process in the ISM. 
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